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Abstract. We propose the use of coherent control of a multi-qubit -cavity QED 
system in order to cxplore novel phase transition phcnomcna in a general class of multi- 
qubit-cavity systems. In addition to atòmic systems, the associated super-radiant 
phase transitions should be observable in a variety of solid-state experimental systems, 
including the technologically important case of interacting quantum dots coupled to 
■ an optical cavity mode. 
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1. Introduction 



There is much current interest in the use of coherent control in order to generate 
novel matter-radiation states in cavity QED and atom-optics systems (J. In addition, 
the field of cavity QED has caught the interest of workers in the field of solid-state 
nanostructures, since effective two-level systems can be fabricated using semiconductor 
quantum dots, orgànic molècules and even naturally-occuring biological systems such as 
the photosynthetic complexes LHI and LHII and in biological imaging setups involving 
FRET (Fluoresence Resonance Energy Transfer) [2]. Such nanostructure systems could 
then be embedded in optical cavities or their equivalent, such as in the gap of a photonic 
band-gap material We refer to Ref. |lj for a discussion of the size and energy-gaps 
of the artificial nanostructure systems which can currently be fabricated experimentally. 

In a parallel development, phase transitions in quantum systems are currently 
attracting much attention within the solid-state, atòmic and quantum information 
communities [3 [üj [TJ [HJ - Most of the focus within the solid-state community has been 
on phase transitions in electrònic systems such as low-dimensional magnets E] while 
in atòmic physics there has been much interest in phase transitions in cold atom gases 
and in àtoms coupled to a cavity. In particular, a second-order phase transition, from 
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normal to súper radiant, is known to arise in the Dicke model which considers N two-state 
àtoms (i.e. 'spins' or 'qubits' pJIE]) coupled to an electromagnètic field (i.e. bosonic 
cavity mode) [0] EU HI] . The Dicke model itself has been studied within the atòmic 
physics community for fifty years, but has recently caught the attention of solid-state 
physicists working on arrays of quantum dots, Josephson junctions, and magnetoplasmas 
[T5] . Its extension to quantum chaos |14j . quantum information J5| and other exactly 
solvable models has also been considered recently [15] . It has also been conjectured that 
superradiance could be used as a mechanism for quantum teleportation [T7j . 

Here we extend our discussion in Ref. [TH] on the exploration of novel phase 
transitions in atom-radiation systems exploiting the current levels of experimental 
expertise in the area of coherent control. The corresponding experimental set-up can 
be a cavity-QED, atom-optics, or nanostructure-optics system, whose energy gaps and 
interactions are tailored to be the required generalization of the well-known Dicke model 
11 . We show that, according to the vàlues of these control parameters, the phase 
transitions be driven to become first-order. 

2. The Model 

The well-known Dicke model from atom-optics ignores interactions between the 
constituent two-level systems or 'spins' [TT]. In atòmic systems where each 'spin' is 
an atom, this is arguably an acceptable approximation if the àtoms are neutral and the 
atom-atom separation d^> a where a is the atòmic diameter. However there are several 
reasons why this approximation is unlikely to be vàlid in typical solid-state systems. 
First, the 'spin' can be represented by any nanostructure (e.g. quantum dot) possessing 
two well-defined energy levels, yet such nanostructures are not typically neutral. Hence 
there will in general be a short-ranged (due to screening) electrostatic interaction 
between neighbouring nanostructures. Second, even if each nanostructure is neutral, 
the typical separation distance d between fabricated and self-organised nanostructures 
is typically the same as the size of the nanostructure itself. Hence neutral systems 
such as excitonic quantum dots will still have a significant interaction between nearest 
neighbors [TÜ] . 

Motivated by the experimental relevance of 'spin-spin' interactions, we introduce 
and analyze a generalised Dicke Hamiltonian which is relevant to current experimental 
setups in both the solid-state and atòmic communities [20] • We show that the presence 
of transverse spin-spin coupling terms, leads to novel first-order phase transitions 
associated with super-radiance in the bosonic cavity field. A technologically important 
example within the solid-state community would be an array of quantum dots coupled 
to an optical mode. This mode could arise from an optical cavity, or a defect mode 
in a photonic band gap material [20]- However we emphasise that the N 'spins' may 
correspond to any two-level system, including superconducting qubits and àtoms [TÏÏ]l2Ti]. 
The bosonic field is then any field to which the corresponding spins couple fHH^O]- Apart 
from the experimental prediction of novel phase transitions, our work also provides an 
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interesting generalisation of the well-known Dicke model. 

The method of solution that we present here is in fact vàlid for a wider class of 
Hamiltonians incorporating spin-spin and spin-boson interactions [21]. We follow the 
method of Wang and Hioe [TTj, whose results also proved to be vàlid for a wider class 
of Dicke Hamiltonians. We focus on the simple example of the Dicke Hamiltonian with 
an additional spin-spin interaction in the y direction. 

N í A e 1 

H = flta + Ç |_^=( a + ct)(af + aj) + -of - Ja] ■ af +1 j (1) 

N í A e 1 

= a]a + £ |^(« + flt X + ^ - Ja í • • ( 2 ) 

Following the discussion above, the experimental spin-spin interactions are likely to 
be short-ranged and hence only nearest-neighbor interactions are included in H . The 
operators in Eqs. 1 and 2 have their usual, Standard meanings. 



3. Results 



To obtain the thermodynamical properties of the system, we first introduce the Glauber 
coherent states \a) of the field [12] where a\a) = a\a), (a\aJ = (a\a*. The coherent 
states are complete, J — |q;}(q;| = 1. In this basis, we may write the canonical partition 
function as: 

Z(JV,r) = j:/— <s|(a|e-^|a>|8> (3) 

As in Ref. [TT|, we adopt the following assumptions: 

(i) a/y/N and a^/yN exist as iV — > oo; 

(ii) limjv^oo lim^oo J2^=o ( '~^ JV ' > can be interchanged 
We then find 



where 



Z(N,T)= í ^ e -^l 2 Tre-^' (4) 

J 7T 



^ / 2ARe(«) x e z Y Y \ 



We first rotate about the y-axis to give 



v I /2ARe(a)\ 2 / e ^ 2 



*=- J £(^wJ nèK + ^)- («) 

We note here that the resulting hamiltonian is of the type of an Ising hamiltonian with 
a transverse field, and it exhibits a divergence in concurrence at its quantum phase 
transition (see, e.g., [Z]). This particular model is instrumental in understanding the 
nature of coherence in quantum systems. Going back to the calculations, we may now 
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diagonalise H' by performing a Jordan-Wigner transformation, passing into momentum- 
space and then performing a Bogoliubov transformation (see, for example, Ref. jü]). We 
then have, in terms of momentum-space fermion operators 7/-, the diagonalised H': 



N 



k=i 



with 



&(a) = 2JJl + {g{a)Y - 2g(a) 



We may then write 



/ 2ÀRe( 



a 



JVN 



N 



where 



k=l 



H k = &(a)(7Í7fc 



(7) 

(8) 
(9) 

(10) 
(11) 



From the transformation, we may associate the spin-up state with an empty orbit on the 
site and a spin-down state with an occupied orbital. Using the commutation relations 
for the 7fc and the fact that 7fc|0) = (see, for example, Ref. pj, we obtain 



Z(N,T) 



d 2 a 



7T 



N 

f[{e-^ + e 

k=l 



}• 



Writing d 2 a = dRe(a)dIm(a), w = Re(a) and integrating out Im(a) we obtain 

1 í c / u , e ~ /3M,2+ S^i{ log [ cosh (§^( a; ))]+ lo s( 2 )} . 



Z(N, T) 



7Ï 



We now let x = w/y/N. Writing £f =1 as £ J 2n dk, yields 

Ipïr 



Z(N,T) 



dx{e~^ 2+I ^y 



where 



and 



1 f 2 * I 

H x ) = 7T dk \\og 
2n Jo 



cosh l^k(x) 



+ log(2) 



(12) 



(13) 



(14) 



(15) 



Ç k (x) = 2J^l + (g(x)) 2 -2g(x) cosk . (16) 

From here on, we omit the log(2) term in I(x) since it only contributes an overall factor 
to Z(N, T). 

Laplace's method now telis us that 

Z(N,T) oc max exp \N[-(3x 2 + I(x)}\ . (17) 

— oo<x<oo l. J 

Denoting [— (3x 2 + I(x)] by Çl(x), we recali that the super-radiant phase corresponds to 
Q(x) having its maximum at a non-zero x [O]. If there is no transverse field, i.e., if 
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J = 0, and the temperature is fixed, then the maximum of Q(x) will split continuously 
into two màxima symmetric about the origin as À 2 increases. Hence the process is a 
continuous phase transition. 

However the case of non-zero J is qualitatively different from J = 0. As a result of 
the frustration induced by the tranverse nearest-neighbour couplings, there are regions 
where the super-radiant phase transition becomes first-order. Hence the system's phase 
transition can be driven to become first-order by suitable adjustment of the nearest- 
neighbour couplings. This phenomenon of first-order phase transitions is revealed by 
considering the functional shape of íl(x), as shown in Fig. ^ 

Figure El shows the value of x that maximises íl(x) at fixed e and two different 
vàlues of J. From the two lines, we can see that the spin-spin coupling actually acts to 
inhibit the phase transition. As we increase J from 0.8 to 1.0 we can see that the value 
to which we have to increase A to induce a phase transition is higher. 

Figure El plots the maximiser of Q(x) with À fixed at a value of 1.3. For small J, 
the local (non-zero) maximum of íl(x) converges to zero as we increase e and the system 
is no longer super-radiant. This is no longer the case if J is increased. In this case, 
Q(x) has a global maximum when e is small; however as e increases, the non-zero local 
màxima becomes dominant and as a result a first-order phase transition occurs. We 
note that the barriers between the wells are infinite in the thermodynamic limit, hence 
we expect that the sub-radiant state is metastable as e increases. This observation also 
suggests the phenomenon of hysteresis, which awaits experimental validation. 

In Fig. 0]we consider the order parameter of the transition, ( s jf-). Following the 
same method as above, we may calculate this to be equivalent to x 2 with an additional 
■L· term that comes from the imaginary part of the coherent states of the radiation field 
|22j . We can see from the figure that as we lower (3 we drive the system first through 
a first order phase transition and then through a continuous phase transition. Thus we 
are able to achieve both a first and second order phase transition by varying the one 
parameter, f3. 

4. Conclusion 

In conclusion, we have shown that the experimentally relevant spin-spin interaction 
in the Dicke model transforms it into an Ising-hamiltonian with a photon-field 
dependent transverse field, which allows for an existence of both first-order and second- 
order phase transitions as parameters vary. Our results highlight the importance 
of spin-spin coupling terms in spin-boson systems and opens up the possibility of 
coherently controlling the competition between the sub-radiant and super-radiant states 
in experimental atom-radiation systems [21] . 

References 

[1] See, for example, A. Olaya-Castro, N. F. Johnson and L Quiroga, Phys. Rev. A 70, 020301 (R) 
(2004); J. Opt. B: Quantum Semiclass. Opt. 6, S730 (2004); Phys. Rev. Lett., in press (2005). 



Exploring phase transitions of a multi- qubit-cavity system 



6 



[2] A.J. Berglund, A.C. Doherty and H. Mabuchi, Phys. Rev. Lett. 89, 068101 (2002); J. Gilmore and 
R.H. McKenzie, cond-mat/0401444 A. Olaya-Castro, C.F. Lee and N.F. Johnson, in preparation 
(2005). 

[3] P.M. Hui and N.F. Johnson, Sòlid State Physics, Vol. 49, edited by H. Ehrenreich and F. Spaepen 

(Acadèmic Press, New York, 1995). 
[4] A. Olaya-Castro and N. F. Johnson, Handbook of Theoretical and Computational Nanotechnology, 

in press (2005); |qüan^ph7 0406133. 
[5] L.P. Kadanoff, Statístical Physics (World Scientific, Singapore, 2000). 
[6] S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 1999). 
[7] T.J. Osborne and M.A. Nielsen, Phys. Rev. A 66, 032110 (2002); R. Somnia, G. Ortiz, H. Barnum, 

E. Knill, and L. Viola, |quant-ph/0403035l 
[8] M.A. Nielsen and I.L. Chuang, Quantum Computation and Quantum Information (Cambridge 

University Press, Cambridge, 2002). 
[9] R.H. Dicke, Phys. Rev. 170, 379 (1954). 
[10] K. Hcpp and E.H. Licb, Ann. Phys. (N.Y.) 76, 360 (1973). 

[11] Y.K. Wang and F.T. Hioe, Phys. Rev. A 7, 831 (1973); F.T. Hioe, Phys. Rev. A 8, 1440 (1973). 
[12] R. Glauber, Phys. Rev. 131, 2766 (1963). 

[13] T. Vorrath and T. Brandes, Phys. Rev. B 68, 035309 (2003); W.A. Al-Saidi and D. Stroud, Phys. 
Rev. B 65, 224512 (2002); X. Zou, K. Pahlke and W. Mathis, |qüant-ph/02010l"ï| S. Raghavan, H. 
Pu, P. Meystre and N.P. Bigelow, |cond-mat/0010140| N. Nayak, A.S. Majumdar and V. Bartzis, 
J. Nonlinear Òptics 24, 319 (2000); T. Brandes, J. Inoue and A. Shimizu, |cond-mat /9908448 
and cond-ma t/9908447 

[14] C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 (2003); Phys. Rev. E 67, 066203 (2003). 
[15] N. Lambert, C. Emary and T. Brandes, Phys. Rev. Lett. 92, 073602 (2004); S. Schneider and G.J. 

Milburn, |qüant-ph/0112"042| G. Ramon, C. Brif and A. Mann, Phys. Rev. A 58, 2506 (1998); 

A. Messikh, Z. Ficek and M.R.B. Wahiddi n, |quant-ph/0303Ï00| 
[16] C. Emary and T. Brandes, |quant-ph /0401029| S. Mancini, P. Tombesi and V.I. Man'ko, 

|quant-ph/9806034| 

[17] y'n. Chen et al., |cond-mat/0502412 



[18] C.F. Lee and N.F. Johnson, Phys. Rev. Lett. 93, 083001 (2004). 

[19] L. Quiroga and N.F. Johnson, Phys. Rev. Lett. 83, 2270 (1999); J.H. Reina, L. Quiroga and N.F. 

Johnson, Phys. Rev. A 62, 012305 (2000). 
[20] E. Hagley et al., Phys. Rev. Lett. 79, 1 (1997); A. Rauschenbeutel et al, Science 288, 2024 (2000); 

A. Imamoglu et al, Phys. Rev. Lett. 83, 4204 (1999); S.M. Dutra, P.L. Knight and H. Moya- 

Cessa, Phys. Rev. A 49, 1993 (1994); Y. Yamamoto and R. Slusher, Physics Today, Junc (1993), 

p. 66; D.K. Young, L. Zhang, D. D. Awschalom and E.L. Hu, Phys. Rev. B 66, 081307 (2002); 

G.S. Solomon, M. Pelton and Y. Yamamoto, Phys. Rev. Lett. 86, 3903 (2001); B. Moller, M.V. 

Artemyev and U. Woggon, Appl. Phys. Lett. 80, 3253 (2002); N. F. Johnson, J. Phys. Condens. 

Matter 7, 965 (1995). 
[21] C.F. Lee, T.C. Jarrett and N.F. Johnson, unpublished. 
[22] G.C. Duncan, Phys. Rev. A 9, 418 (1974). 



Exploring phase transitions of a multi- qubit-cavity system 



7 



Figure 1. Demonstration as to what the function Q(x) looks like across a first-order 
phasc transition as A and x are varied. Here J = 1.0, e = 1.1 and (3 = 100.0. 




Figure 2. The value of x at which there is a maximum in Çl(x) as A increases, for 
J = 1.0 (dashed line) and J = 0.8 (sòlid line). In both cases e = 1.1 and (3 = 100. 
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Figure 3. The maximiser of Í7 shown as a function of J and e. Here A = 1.3 and 
[3 = 100. 
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Figure 4. Plot of the order parameter, 9 = for the phase transition with 

A = 0.9, J = 1.0 and e = 1.1. 






